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Abstract 

We present detailed calculations of the xt\ coannihilation channels that have the largest 
impact on the relic x density in the constrained minimal supersymmetric extension of the 
Standard Model (CMSSM), in which scalar masses m o, gaugino masses mj /2 and the trilinear 
soft supersymmetry-breaking parameters A 0 are each assumed to be universal at some input 
grand unification scale. The most important t\t* and tfti annihilation channels are also 
calculated, as well as til coannihilation channels. We illustrate the importance of these new 
coannihilation calculations when Ao is relatively large. While they do not increase the range 
of 777-1/2 and hence m x allowed by cosmology, these coannihilation channels do open up new 
‘tails’ of parameter space extending to larger values of mo- 
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1 Introduction 


A favoured candidate for cold dark matter is the lightest supersymmetric particle (LSP), 
which is generally thought to be the lightest neutralino x [HI i n the minimal supersym¬ 
metric extension of the Standard Model (MSSM). It is common to focus attention on the 
constrained MSSM (CMSSM), in which all the soft supersymmetry-breaking scalar masses 
m 0 are required to be equal at an input superysmmetric GUT scale, as are the gaugino 
masses mm and the trilinear soft supersymmetry-breaking parameters Aq. These assump¬ 
tions yield well-defined relations between the various sparticle masses, and correspondingly 
more de fini te predictions for the relic abundance fl~h 2 and observable signatures. This paper 
is devoted to relic-abundance calculations including coannihilations of the lightest neutralino 
X with ti , the lighter supersymmetric partner of the top quark ||. 

The range 0.1 < fl~h 2 < 0.3 is generally thought to be preferred by astrophysics and 
cosmology ||. Lower values of hi-fi 2 might be possible if there is some other source of 
cold dark matter, but higher values are incompatible with observation. The regions of the 
mi/ 2 , mo plane where the relic density falls within the preferred range 0.1 < fl~fi 2 < 0.3 have 
generally been divided into four generic parts. There is a ‘bulk’ region at moderate mm 
and mo JI]]. Then, extending to larger m. 1 / 2 , there is a ‘tail’ of the parameter space where 
the LSP x is almost degenerate with the next-to-lightest supersymmetric particle (NLSP), 
which is in this region the fi, the lighter supersymmetric partner of the r lepton. Along 
this ‘tail’, efficient coannihilations 0, i! keep fl~fi 2 down in the preferred range, even for 
larger values of m~ S 1,0- 0- At larger mo, close to the boundary where electroweak 
symmetry breaking is no longer possible, there is the ‘focus-point’ region where the LSP has 
a larger Higgsino component and m~ is small enough for fl~fi 2 to be acceptable [|ll|]. Finally, 
extending to larger mm and m 0 at intermediate values of m^/mo, there may be a ‘funnel’ 
of CMSSM parameter space where rapid direct-channel annihilations via the poles of the 


heavier Higgs bosons A and H keep hi-fi 2 in the preferred range [|1|, [13 


I 11 this paper, we emphasize the significance of coannihilation of the LSP x with U, the 
lighter supersymmetric partner of the t quark 0. This mechanism opens up another ‘tail’ 
of parameter space, this time extending to larger values of mo- It is not relevant for the 
small values of A 0 considered in previous coannihilation calculations E, 0- 0. but may 
be important for large A 0 , as we demonstrate in this paper. Coannihilations of x with t\ 
are important when the latter is the NLSP, just as x.f\ coannihilations are important when 
the t 1 is the NLSP. In the latter case, one must also consider coannihilations with the e\ 
and fli, which are not much heavier than the f\ i 1- I- ITOj . There are also regions of 
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CMSSM parameter space where both the t\ and f\ are close in mass to the LSP x, and 
t\T\ coannihilations must also be considered. We present here detailed calculations of the 
matrix elements and cross sections for all the leading xti and t\i coannihilation processes, 
and illustrate their importance for f2~h 2 in some instances in the CMSSM when A 0 ^ 0. 

The structure of the paper is as follows. In Section 2 we recall some important features 
of LSP relic-density calculations in general, and coannihilations in particular. Then, in 
Section 3 we compare the relative magnitudes of the xX-> Mi*'* and t\fS*) processes 
for some specific choices of the CMSSM parameters. Section 4 provides an overview of 
the implications of xti coannihilation and related processes for the regions of the m\/ 2 ,m o 
plane allowed by the constraint 0.1 < Vt-ti 2 < 0.3 for various choices of the other CMSSM 
parameters. Relevant details of our calculations of the matrix elements are contained in an 
Appendix. 


2 Formalism for Annihilation and Coannihilation 


The density of neutralino relics left over from the early Universe may be determined relatively 
simply in terms of relevant annihilation cross sections, using the Boltzmann rate equation to 
determine a freeze-out density. The relic density subsequently scales with the inverse of the 
comoving volume, and hence with the entropy density. In the MSSM framework discussed 
here, since neutralinos are Majorana fermions, the S'-wave annihilation cross sections into 
fermion-antifermion pairs are suppressed by the masses of the final-state fermions, and it is 
therefore necessary to compute P -wave contributions to the cross sections |l|. 

The rate equation for a stable particle with density n is 


^ - (fTUrei) (n 2 - n 2 eq ) , (1) 

where n eq is the equilibrium number density and (<7n re i) is the thermally averaged product of 
the annihilation cross section a and the relative velocity v rfi \. In the early Universe, we can 
write R/R = (87tGatp/ 3) 1//2 , where p = 7r 2 g(T)T 4 /30 is the energy density in radiation and 
g(T) is the number of relativistic degrees of freedom. Conservation of the entropy density 
s = 27r 2 h(T)T 4 /45 implies that R/R = —T/T — h'T/3h where h' = dh/dT. Generally, we 
have h{T ) rs g(T). Defining x = T/m and q = n/T 3 h, we can write 


J/. = 171 (tU 1 G n9) 1/2 (h + \rnxh!') (av Tel )(q 2 - q 2 eq ) . (2) 


The effect of the b! term was discussed in detail in 


14], and is most important when the 


mass in is between 2 and 10 GeV. Since we only consider neutralinos that are significantly 
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more massive, we neglect it below (though it is not neglected in our calculations). In the 
case of the MSSM, freeze-out occurs when x ~ 1/20, and the final relic density is determined 
by integrating ([!]) down to x — 0, and is given by 

p~ = mg(0)/r(0)T 0 3 . (3) 

When coannihilations are important, there are several relevant particle species i, each 
with different mass, number density ni and equilibrium number density n eqi . Even in such a 
situation ||, the rate equation ([I]) still applies, provided n is interpreted as the total number 
density, 

n = J2 n i ^ ( 4 ) 

i 

n eq as the total equilibrium number density, 


n, 


eq 


= E 




eq,z 5 


(5) 


and the effective annihilation cross section as 

n, 


(^efWrel) — ^ ) 


eq,z^eqj 


nt 


rel) • 


( 6 ) 

ij ■ “eq 

In (|2j) , m is now understood as the mass of the lightest particle under consideration. For 
T < 77ii, the equilibrium number density of each species is given by |TI], [15 


n, 


eq,* 


= 9i 


d 3 p 


(27,)' 


E/T 


girrijT 


= <h 


2n 2 


K 2 {m i /T) , 


exp (-rrii/T) 


15 T 

1 + tor + '" 


(7) 


where g t is a spin and color degeneracy factor and K 2 (x) is a modihed Bessel function. 
We make the approximation of Boltzmann statistics for the annihilating particles, which is 
excellent in practice. 

We now recall how to compute (ai 2 U re i} for the process 1+2 —> 3+4 in an efficient manner. 
Suppose we have determined the squared transition matrix element \T\ 2 (summed over final 
spins and averaged over initial spins) and expressed it as a function of the Mandelstam 
variables s, t, u. We then express \T\ 2 in terms of s and the scattering angle d CM in the 


center-of-mass frame, as described in 


w(s 


1 

4 


d 3 p 3 


J. We now define 

d 3 p4 


(2n) 3 E 3 (2? t) 3 E 4 

Pais) f +1 


(2vr) 5 (pi +P2-P3 ~Pa) \T\ 


32t r s 1 / 2 


dcos^cM | T 


/-i 


( 8 ) 
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In terms of w (s), the total annihilation cross section a 12 (s) is given by a 12 (s) — w(s)/s 1/l2 pi (s) []. 

The above analysis is exact. To reproduce the usual partial wave expansion, we expand 
|T | 2 in powers of p\{s)/m\. The odd powers vanish upon integration over d CM , while the 
zeroth- and second-order terms correspond to the usual S and P waves, respectively. Each 
factor of pi(s) is accompanied by a factor of cos 9 CM , so we have 

£bco S « OM T| 2 = (|T| 2 „ fcM ^ +1/vI + mL, ou ^_ 1/V5 ) + O(pi) . (9) 


We can therefore evaluate the S- and P -wave contributions to w(s) simply by evaluating 
|T | 2 at two different values of cos# CM ; no integrations are required. 

The proper procedure for thermal averaging has been discussed in [ p~4| , [L5|] for the case of 
m i = and in |TB], Bj for the case of rri\ 7 ^ m 2 , so we do not discuss it in detail here. One 
finds 

(of 2 iYei) = a l2 + b 12 x + 0(x 2 ) , (10) 


where x = T/mi (assuming m\ < m 2 ). In our case, we extract Gq 2 and b\ 2 from the transition 
amplitudes listed in the Appendix for each final state. We set x = 0 to get a 12 , and then 
compute b \2 by setting x to a numerical value small enough to render the 0(x 2 ) terms 
negligible. We then compute a e s and b e g by performing the sum over initial states as in (| 6 |), 
and integrate the rate equation (|2|) numerically to obtain the relic LSP abundance. To a fair 
approximation, the relic density can simply be written as [llj, |J 


Qh 2 


10 9 GeV ” 1 

1 To "> 

g/ M p i(a eS + b eS x f /2)x f 


( 11 ) 


where the freeze-out temperature 2 } ~ m~/ 20 , and g f is the number of relativistic degrees 
of freedom at 2 }. 

This implies that the ratio of relic densities computed with and without coannihilations 
is approximately 


R = 


Q 



( 12 ) 


where a = a + bx /2 and sub- and superscripts 0 denote quantities computed ignoring coanni¬ 
hilations. The ratio x 9 / x s ~ l+^° ln(, 9 e frcr e ff/rq(To), where g e g = J2i gi(mi/mi) 3 ^ 2 e~^ mi ~ mi ^ T . 
For the case where the t\ and \ are almost degenerate, g e g ~ J2i9i — 8 an d x®/x f ~ 1.2. 

1 Our w(s) is also the same as w(s) in ill which is written as W/A in |(J. 
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3 Coannihilation Rates for t\ in the MSSM 


We now use the above formalism to estimate the relative importance of the t\X coannihilation 
processes, t\t\ and t\t\ annihilations calculated in the Appendix. We also take into account 
the yl coannihilations calculated previously 01 and, for completeness, include the t\£ and 
t\l* coannihilations also calculated in the Appendix. 

To compute the effective annihilation cross sections for light sparticles in the MSSM, we 
allow the index i in (|) to run over ti, t\, A, t*, e R , e*, /i I{ and /f*, as well as y. The following 
is the change in a e s compared with || , where 49 of the <ly,- in (|J) were already included: 

Aa eff = 2 (a hil + a^~*)r| +4a~ h r~r h + 8 (<r ilSR + a t - lg . )r t - 1 rg H + 4 (a hfl + a hf ,)r~ r~ tl (13) 

where ?y = n eqi i/n e q . We have taken the e R and /i R (but not the f) to be degenerate in mass, 
thus accounting for the 81 possible initial state combinations. Note that we have summed 
over color states in the cross sections amplitudes listed in the Appendix, and we have taken 
the stop degeneracy factor gi = 3. We list in Table [I] the sets of initial and final states for 
which we compute the annihilation cross sections, using the transition amplitudes given in 
the Appendix. We use q to denote the four light quarks, which we have taken to be massless. 

Table 1: Initial and Final States for t\ Annihilation and Coannihilation Processes 


Initial State 

Final States 


99,19, Zg , **, bb, qq , gh, gH , Zh, ZH , ZA, W ± H^ } 


hh, hH , HH , AA, hA, HA , H+H~ 

* 1*1 

tt 

X*i 

tg, tZ, bW + , tH , th , tA, bH + 

tj 

U, bv 

tj* 

tl 


In the CMSSM, the diagonal entries of the squark mass matrix tend to pick up large 
contributions from the gaugino masses, m 2 3 LLRR 3 0(6)rriy 2 , thus making the squarks heavier 
than the neutralinos. The off-diagonal entry for an up-type squark 0 

m LR = ~m q (A q + /icot/3) (14) 

can, however, be large, particularly for the stops, or for sbottoms at large tan/3 £]. When 
A t is sufficiently large, the lighter stop, t±, can become degenerate with (or lighter than) 

2 Note here our sign convention for A q . 

3 For down-type squarks, the factor cot (3 in (0) is replaced by tan (3. 
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Figure 1: The separate contributions to the xt\ coannihilation cross sections a = a + \bx 
for x = T/m~ = 1/23 as functions of m 0 for (a) my 2 = 230 GeV, T 0 = 1000 GeV and (b) 
my 2 = 450 GeV, A 0 = 2000 GeV. Also shown are the total cross section and , for comparison, 
the much smaller total cross section for xx annihilation. 


the neutralino. Thus, it is when To is large that we expect x.t\ coannihilations to become 
important. 

It is important to distinguish between the effective low-energy parameters A t , etc., 
and the high-energy input parameter A 0 , which are related through the running of the 
renormalization-group equations. For example, for tan/3 = 10 and m 0 = 300 GeV, t\X 
coannihilations are important when my 2 = 200,450,and 670 GeV and To = 1000,2000, 
and 3000 GeV, but these values correspond to A t — 565,1200, and 1700 GeV respectively. 
Furthermore, the values of T for the light squarks are different and typically larger than T t . 

We display in Fig. |1| numerical values of the contributions to a = a + bx/2 (see ©) 
in xM coannihilation, for the representative values x = 1/23, tan /3 — 10, p > 0 and (a) 
mi /2 = 230 GeV, T 0 = 1000 GeV and (b) my 2 = 450 GeV, T 0 = 2000 GeV as functions of m 0 . 
For comparison, the total cross section for XX annihilation to all final states is also shown, as 
a thick dotted line. We see that the xt\ —*■ tg and th coannihilation cross sections dominate 
by large factors over the total XX annihilation cross section, suggesting that they may have a 
greater importance than that suggested by simply comparing Boltzmann suppression factors. 
The feature in Fig. |l](a) at m 0 ~ 400 GeV is due to the threshold for the production of th 
final states. At smaller values of m 0 , this final state is kinematically forbidden. 

Fig. |2| displays similar plots for tyt* annihilation, for the same parameter choices as in 
Fig-B In this case, the dominant tyt\ annihilation cross sections are into gg and hh, and even 
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Figure 2: The separate contributions to the tit* annihilation cross sections cr = a + \bx 
for x = T/m~ = 1/23, as functions of mo for (a) = 230 GeV, Aq = 1000 GeV and (b) 

mi /2 = 450 GeV, A 0 = 2000 GeV. Also shown are the total cross section and, for comparison, 
the much smaller total cross section for xx annihilation. 


subdominant cross sections such as 7 g, Zg, gh and the various quark-antiquark channels are 
far larger than the total xx annihilation cross section. Once again, when Ao = 1000 GeV 
we see thresholds, in this case corresponding to hh production at mo ~ 180 GeV and tt 
production at m 0 ~ 330 GeV. 

The t\ti annihilation cross sections shown in Fig. show that the cross section for anni¬ 
hilation into the tt final state, when it is kinematically open, is also far larger than the total 
XX annihilation cross section. 

A complete study of coannihilation effects must include not only the xf processes consid¬ 
ered previously 0 . 0 . and the xU processes considered above, but also £ 1/1 coannihilations. 
Accordingly, the final set of coannihilation cross sections we present are those for GG and 
t[ti, shown in Fig. f|. We see that, when (a) A 0 = 1000 GeV, the tr, tf and bv e final states 
are the most important, followed by te, bu T and te, whereas (b) the tf and te final states 
are relatively much less important when Ao = 2000 GeV. In all panels of Fig. 0. there are 
coannihilation cross sections much greater than the total xx annihilation cross section, which 
is also plotted. 

The basic reason for the relatively small magnitude of the XX annihilation cross section is 
that it is dominated by the P- wave suppressed cross sections for XX annihilation to fermion 
pairs. This was also the basic reason why xl coannihilation processes were previously found 
to be so important 0, | 8 |, [9], [L(J . 


7 













> 

Qi 

o 


X 

,T 

I'M 

I'M 

<b 




m 0 (GeV) m 0 (GeV) 


t t 


°3CX->all 


Figure 3: The i\i\ —> tt annihilation cross sections a = a + \bx for x = T/m~ = 1/23, 
as functions of m 0 for (a) my 2 = 230 GeV, Aq = 1000 GeV and (b) my 2 = 450 GeV, 
Ao = 2000 GeV. Also shown, for comparison, is the much smaller total cross section for 
XX annihilation. 
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Figure 4: The separate contributions to the l\t\ coannihilation cross sections a = a + |bx 
for x = T/m~ = l/23 ; as functions of for (a) my 2 = 230 GeV ; Ao = 1000 GeV and (b) 
m i /2 — 450 GeV, Aq = 2000 GeV. Also shown, for comparison, is the much smaller total 
cross section for XX annihilation. 















The contributions of the various annihilation channels to cr e g are weighted by the rela¬ 
tive abundances of the y, ti and l\. For a stop degenerate with the y, tit\ annihilation, 
t\ti annihilation and yfi coannihilation are clearly the dominant contributions to Aa e ff, and 
hence to cr e fj in (jBj), and the final neutralino relic density is greatly reduced. As the stops 
become heavier than the neutralinos, their number densities are exponentially suppressed 
and the stop contributions to a e g become less important. Fig. || shows the sizes of the 
separate contributions to a e s from yy annihilation, t \y coannihilation and t\t\, t\t\ anni¬ 
hilations (combined), as functions of the mass difference between the t\ and y. In Fig. || 
we have fixed m 0 = 300 GeV, tan /? = 10, /r > 0, A 0 = (a) 1000 and (b) 2000 GeV, and 
computed a e g for varying which amounts to varying the fractional mass difference 

AM = (rrifj — m~)/m~. For these choices, the stau mass, m^, is much larger than m The 
thin dotted lined is the value of a that one would compute if one ignored all coannihilation 
contributions []. Note that, in the case of close degeneracy between the y and t \, the t\ti and 
tit\ annihilations dominate <f e fj. However, since these contributions are suppressed by two 
powers of n eq ^, they drop rapidly with AM, and neutralino-stop coannihilation takes over. 
For Aq = 2000 GeV, this occurs at AM > 0.18. This contribution in turn falls with one 
power of n eq q, and yy annihilation re-emerges as the dominant reaction for AM > 0.25. 
When AM > 0.35, t\ coannihilation effects can be neglected. In We see the presence of 
kinematic thresholds also in Fig. [5]. In panel (a), we see the threshold for producing a single 
top quark in ty coannihilations, and in panel (b) we see the threshold for producing a ti pair 
in yy annihilations. 

4 Implications of l\ Coannihilations for the Region of 
CMSSM Parameter Space Favoured by Cosmology 

We now explore the consequences of H coannihilations for the region of CMSSM parameter 
space in which 0.1 < h 2 < 0.3, as favoured by cosmology. We display in Fig. |] the 
( m i/ 2 , m o) planes for tan/? = 10 and /i > 0, for the different values of A 0 = (a) 0, (b) 1000, 
(c) 2000 and (d) 3000 GeV. The very dark (red) shaded regions have rrif i or rnq < m~, 
and hence are excluded by the very stringent bounds on charged dark matter [If]. The light 
(turquoise) shaded regions correspond to the preferred relic-density range 0.1 < h 2 < 0.3. 
The dark (green) shaded regions are those excluded by measurements of b — > sy. The 
intermediate (pink) shaded regions in panels (a) and (b) are those favoured by the BNL 
4 This differs from <r~~ because of the number-density weighting factor. 
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Figure 5: The separate contributions to the total effective cross sections <j e ff for x = T/m~ 
1/23, as functions of AM = (rrif H — obtained by varying with (a) H 0 

1000 GeV and (b) A 0 — 2000 GeV, both for tan (3 = 10, p > 0 and m 0 = 300 GeV. 


measurement of g lt — 2 at the 2-a level []T7|. The (black) dashed line in panel (a) is the 
contour m~± = 104 GeV, corresponding to the kinematic range of LEP 2, and the (red) 
dotted line is the contour mh = 114 GeV, as evaluated using FeynHiggs |18|, corresponding 
to the LEP lower limit on the mass of the Higgs boson. These contours are also glimpsed in 
the other panels. 

Panel (a) of Fig. is indistinguishable from analogous plots shown previously 
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ti coannihilations are not important when Ho = 0. Panel (b) shows how t\ coannihilation 
generates a ‘tail’ of allowed CMSSM parameter space, extending the ‘bulk’ region of the 
(mi/2,m 0 ) plane up as far as m 0 ~ 450 GeV. The boundary of the region where < m~, 
and hence the tfx coannihilation region, slopes to the left as m 0 is increased, since we have 
fixed Ho over the plane. As mo is increased, the impact of the H-dependent off-diagonal term 
is diminished and the light stop is heavier. To compensate for this, one must decrease m \/2 
to obtain the necessary degree of degeneracy between t\ and y. This effect is also seen in 
panels (c) and (d) of Fig. []. 

In the particular cae of Ho = 1000 GeV shown in panel (b) of Fig. the t\ coannihilation 
tail happens to be excluded by b —> sy. However, this is not the case for (c) H 0 = 2000 GeV 
and (d) H 0 = 3000 GeV, where the t) coannihilation tail extends up to m 0 ~ 900 and beyond 
1350 GeV, respectively. For (b) H 0 = 2000 GeV, the region favoured by — 2 is hidden 
on the left, inside the excluded region. However, around H 0 = 1500 GeV there is a xt\ 
coannihilation region that satisfies both the b —> sy and — 2 constraints. 

The cosmologically-favoured region is broadened significantly when tt production is kine- 
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Figure 6: The (mi/2,m 0 ) planes for tan (3 — 10, /jl > 0 and A 0 = (a) 0, (b) 1000, (c) 2000 
and (d) 3000 GeV. The very dark (red) shaded regions are excluded because the LSP is the 
t\ or the f]. The dark shaded regions are excluded by b —> s'y. The intermediate shaded 
regions are favoured by g /t — 2 at the 2-a level. The light (turquoise) shaded regions are those 
favoured by cosmology, with 0.1 < h 2 < 0.3 after the inclusion of coannihilation effects. 
The (black) dashed lines are the contours m~± = 104 GeV, and the (red) dotted lines are the 
contours mh = 114 GeV. 
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matically allowed in yy annihilation The t-channel stop exchange contribution to this process 
is significantly enhanced when t\ is relatively light: my ~ m~. This feature can be seen in 


Fig. 0(c) and Fig. 0(d). The t\t coannihilations are only important in the corner area near 
the instep of the dark (red) shaded region where rrif x ~ my ~ m~. 

Fig. 0 shows analogous (mi/ 2 ,m 0 ) planes for the larger values tan (3 = 20 and 30, both 
again for y > 0. We see in panel (a) that the t\ coannihilation tail extends up to m o ~ 
900 GeV when tan/3 = 20 and A 0 = 2000 GeV, beyond the region forbidden by b — > sy. 
We then see in panel (b) that the ti coannihilation tail extends beyond m 0 ~ 1400 GeV 
when tan/3 = 20 and Aq = 3000 GeV. A similar portion of the t\ coannihilation tail that 
is consistent with b — > sy is also visible in panel (c), for tan/3 = 20 and Aq = 2000 GeV, 
where mo ~ 900 GeV. For fixed A 0 , as tan/3 is increased the b —► sy constraint becomes 
more severe. 

The t\ coannihilation tails do not increase the allowed range of m.1/2, as was the case for 
the l coannihilation tail E11 - 0 , but they do add a significant filament to the CMSSM 
region allowed by experiment and favoured by cosmology. 

I11 the above illustrative examples, we have chosen the sign of Aq to be the same as that 
of n, with the effect of maximizing the stop off-diagonal mass terms, and hence minimizing 
my. If the sign of A 0 is switched while keeping /i > 0, an analogous xti coannihilation 
region is found only at even larger negative A 0 . However, the constraint is more severe 
for A 0 < 0, and also for /j, < 0, excluding the region of parameter space of interest in the xh 
coannihilation context. 

We have not discussed in this paper the potential constraint imposed by the absence 
of colour and charge breaking (CCB) vacua, or at least the suppression of transitions to 
CCB vacua. This constraint would restrict the value of Aq relative to mo |19|. If the CCB 


constraint is imposed, some of the parameter space for small mo will be excluded, but there 
will still be regions at high m 0 where xt\ coannihilation is crucial. The issue of compatibity 
with the potential g IJr — 2 constraint would then arise. Since raising mo suppresses the 
neutralino-proton elastic scatttering cross section, xt\ coannihilation would tend to lower 
the neutralino direct detection rate for the same range of m~. 


5 Conclusions and Open Issues 

We have documented in this paper the potential importance of xti coannihilation in delin¬ 
eating the preferred domain of CMSSM parameter space for A 0 7^ 0. In this paper, we have 
only sctratched the surface of this subject, whose higher-dimensional parameter space merits 


12 










m 0 (GeV) m 0 (GeV) 


tan B = 20, An = 2000 GeV, u > 0 tan B = 20, An = 3000 GeV, u > 0 



200 300 400 500 600 700 800 300 400 500 600 700 800 


m l/2 (GeV) 

m l/2 (GeV) 

tan B = 30, An = 2000 GeV, u > 0 tan B = 30, An = 3000 GeV, u > 0 



m l/2 (GeV) 

m l/2 (GeV) 


Figure 7: Notation as for Fig. |/| Panels (a,b) are for tan (3 = 20, /x > 0, and panels 
(c,d) are for tan/? = 30,/r > 0. Panels (a,c) are for A 0 = 2000 GeV, and panels (b,d) for 
4, = 3000 GeV. 
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more detailed exploration. The Appendix provides details of the diagrammatic calculations 
that should be sufficient for our results to be verified and used by other authors. Although 
applied in the context of the CMSSM, our results may also be used in more general MSSM 
contexts. However, other coannihilation processes are also important in other regions of the 
general MSSM parameter space. For example, in the CMSSM the sbottom mass is generally 
larger than the stop mass even for large tan (3. However, if one allow non-degeneracy in 
the scalar soft breaking mass term, a sbottom NLSP becomes possible m • A complete 
calculation of the LSP relic density in the MSSM requires a careful discussion of all such 
coannihilation possibilities. 
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Appendix 

We calculate the tree contributions to the coannihilation and annihilation amplitudes that 
are leading order in a s , a t = A 2 / 47 t and/or = A 2 / 47 t. We keep terms that are of the same 
order in y/a s a em: A Joi t oi em , a.\,a. em , ^a s aw, and/or ^a^aw, but neglect most terms 

in a em and/or aw- This approximation is motivated by the relatively large magnitudes of 
a s and a t , and the fact that a& ~ a t at large tan/?. We have however included s-channel 
Z, and 7 exchange for tit* —> tt and checked that they make negligible contribution to the 
total cross section and therefore neglect these channels in the remaining amplitudes. 

Couplings 

Here we list couplings used in the calculation. Factors not written, such as i, y’s, QCD 
generators and momentum , are taken into account in the calculation of amplitudes squared, 
which are listed subsequently. 
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Squared Amplitudes 

Below is the list of the amplitudes squared. Note that, for identical-particle final states, one 
needs to divide them by two when performing the momentum integrations. 
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27n 1 x7yn + 2/ 4 ^ TiyXTy + 27iyX7yi + 27iyX7yn + 27yX7yi + 

i=l 

27yx7yn + 27yix7yn (A4) 
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t\t\ —> bb 


i- channel 

h annihilation 


s-channel H annihilation 


. t-channel Xn 2 ) exchange 


, s-channel gluon annihilation 


h 

— Cl 1 _ tl _ h Cb-b-h 


f'2 

= Cb-b-H 


fzLL 

r'iL s~iL 

— Or l ~+ Or j ~+ 

tl-b-xT tl-b-xT 


fl,LR 

r'iL r^R 

— O 7 , ~+ Lx 7 , ~+ 

ti-b-xl ti-b-xj 


f’iRL 

s~iR s~iL 

— Or j ~+ U: , — 

ti-b-x7 ti-b-xj 


firm 

r'lR r^R, 

— Or , ~+ U: , ~ + 

ti~b~xj *1 -b-xj 


h 

= Cf 1 _l 1 _ g Cb-b-g 



= 2(s-4m 2 )/{s-rn 2 h ) 2 


TnxT n 

= 2(s-4m 2 b )/(s-m 2 H ) 2 


T m xT m 

= ((flmifiRR + /3Ll/3Ll)( t 

{2m 2 - 


+/slr/3Rl) m^+ { s - 

2m^) - 


IIrrIIrl) ™b {m 2 h - m 2 b - t) 


flRLfluf) m xj m b Of, - ~ t) - 


(flrmflrm + flm.flmiy 

< m xt 

TjvXTjy 

= (2/9)(2s(s-4mf 1 )-2(w- 

- tf)/s 

TixTu 

= 2(s — 4m^)/((s — m 2 H )(s - 

- ml )) 

TixT m 

= ((fin + /W.)2m fc (t + 

-™l) 





= 0 


TnxTui 

= (RIll + flRR)2m b (t + ml 



( m l 


m 


b t)~) + (flRLfzLR 


l 2 .+ 

xt 


+ (fsRL + /3Li?) m v+( S ^ m b)) 


(«- Am D) 


/((s-m 2 H )(t-m 2 -+)) 


T u xT w = 0 
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T m xT lv 


0 


2 

\T \ 2 = fi%xTi + f 2 T u xT u + E TuixTui + f^TiyxTiy + 2/i/ 2 7fx7)i + 

i,j =1 

2 2 

2/i ^ 7fx7in + 2 /i/ 4 7ix7iv + 2/ 2 7nx7m + 2j 2 j i {Ti\xT\^ + 

2=1 2=1 

2 

2/ 4 E 7 ni> < ^v (A5) 

2=1 

tit/ —> qq, q = u, d, s, c 

I. s-channel gluon annihilation 

Cq-q-g 

(2/9) (2 s(s — 4m|) — 2 (w — t) 2 )/s 2 
fi%x% (A6) 

Ml —* M 

I. t-channel ti exchange 

II. u-channel / 4 exchange 


/l — Cti-tl-9 [Ch-ti-s 

A = Cfi-h-g Qi-ti-J? [Qi-ti-fl 


71x71 = 

(-4/9)(2(m?+t))/(t-m?) 2 


7iix7ii = 

(-4/9)(2(m|+u))/(u-m|) 2 


7ix7) t = 

(-4/9) (6m? - u - t - 2s)/((t - m|)(u - m?)) 


T 2 = 

7/71x71 + flTuxTu + 2f 1 f 2 T 1 xT u 

(A7) 


t~it~/ —» Z/i [Ztf] 

I. t-channel t( 1)2 ) exchange 

II. u-channel t( 1)2 ) exchange 


A = 

7)x7I = 
ITI 2 = 
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fl — Qi -ti-Z Ck-U-h 

ri _ s^i 

J 2 — Wi -ti-Z ° h-U-h 

Ti*Ti = - ( 2m i + 2t ~ m l - (ml - t) 2 /m|) /((t - m?) (t - m?)) 

^ix^i = “(2m? + 2u - m| - (m? - u) 2 /m 2 z )/{{u - m\)(u - m? )) 

= - (4m? — rn\ — s — (m| - t) (u - m? ) /m|) / ((t - m? ) (« - m? )) 


2 2 



i=l i—1 




I. s-channel h annihilation 

II. s-channel H annihilation 

III. t-channel t 2 exchange 

IV. u-channel t 2 exchange 


fi 

/2 

h 

u 


%x% 

T n xT u 

T m xr m 

TiyXTiv 

TixTn 

%xTiv 

r n xT m 

TuxTiv 


C^-h-h Ch-Z-A 
Cir-h-H C'H-Z-A 


Qi- 

—t2~Z —t^—A 



C'fi- 

—t2~Z ^t\—t2—A 



((* 

2 \2 / 2 
~ m A) ! m Z ~ 

(2s- 

m| + 2 m 2 A ))/(s - ml) 2 

((a 

2 \2 / 2 
~ m A) t m Z ~ 

(2s - 

m| + 2 m\))/(s - m 2 H ) 2 

(( mi-t) 2 /m 2 z - 

( 2m l 

+ 2t))/(t-ml) 2 


((ml - u) 2 /m 2 z - (2 ml + 2 u))/(u - m|) 2 

((s - m 2 A ) 2 /m 2 z - (2s - m| + 2m^))/((s - m? H )(s - m£)) 

((s - m\)(ml - t)/m 2 z - m| + 2m^ - «)/((s - m 2 h ){t - m|)) 
((s - m\)(ml - u)/m 2 z - m\ + 2 m\ - t)f((s - rn 2 h )(u - m|)) 
((s - m^)(m? - t)/m| - m? + 2m^ - w)/((s - m^)(t - m|)) 
((s - m\)(ml - u)/m 2 z - m? + 2m^ - t)/((s - m%)(u - m|)) 


(A8) 
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T m xT w 


= ((u ~ m 2 h ) {ml - t )/m| - (4m? - m\ - s))/{(t - m|) (u - m|)) 

|T| 2 = rtT 1 xT 1 + flT n xT u + flT m xT m + fjT w xT w + 2f 1 f 2 T l xT ll + 
2/i/3^ix7ih + 2fif 4 TixTiy + 2f 2 f 3 TuxTm + 2f 2 j\T n x%y + 

S/s/^iixTiv (A9) 

[f+w-] 

I. s-channel h annihilation 

II. s-channel H annihilation 

III. t-channel &(i, 2 ) exchange 


fi 

f'2 

n 

TiX% 


TwxTm 

TjxT u 

TjxTm 

^nX^m 

ITI 2 


Qi-h-h C h -W+-H- 
C H -W+-H- 
C'ti-bi-W+ Ct 1 -bi-H+ 

(( S- - m 2 H+ ) 2 /m 2 w - {2s - m 2 w + 2m 2 H +))/{s - m 2 h ) 2 

((s - m 2 H+ ) 2 /m 2 w - {2s - m 2 w + 2m 2 H+ ))/{s - m 2 H ) 2 

(« - - 2 {m 2 k + *))/((* - - m ? )) 

((s - m 2 H+ ) 2 /m 2 w - {2s - m^ + 2m|+))/((s - m^)(s - m 2 h )) 

((s - m|+)(m| i - t)/m^ - (m| + 2m|+ - «))/((s - m 2 h ){t - mg.)) 

((«■ - m H+)( m i 1 - - (' m l + 2-m^+ -u))/{{s- m 2 H ){t- mg.)) 

2 2 

/i M + f 2 T u xT u + EE fitiTwxTm + 2f 1 f 2 T 1 xT 11 + 

i=l jr = l 

2 £ j\fl%xT m + 2 £ f 2 flT u xT m (A10) 


—> HH [hh] [hH] 


I. t-channel t(i ; 2 ) exchange 

II. n-cha nn el t( 12 ) exchange 

III. point interaction 

IV. s-channel h annihilation 



V. s-channel H annihilation 


n 

fi 

fa 

fa 

fa 


faxT, 

T n xTu 

TmxT m 

fayxfay 

TyXTy 

faxTn 

TixTm 

faxT lw 

faxTy 

r n xr m 

T u xT iy 

T a xTy 

T m xT w 

T m xTy 

fayXTy 

\T\ 2 


[Cf^ti-h] [ c h-ti-h Qi-ti-n] 
Ch-ii-H-H [^h-ti-h-h] fa'h-h-h-H] 


- 

° ti-U-H 

r? - 

Ly ti-ti-H 


Ch-ti-h Ch-H-H [Qi-h-Zi Ch-h-h] [C^-ix-h Ch-h-H. 


Cti-h-H C H-H-H [Cil-h-H Ch-h-H] [Cti-h-H Ch-H-H. 


l/((u-mf.)(u-m|.)) 

1 

l/(a - <? 

1 fas - mj,) 2 


l/((t - ml)(s - m 2 h )) 

l/({t-m\)(s-m 2 H )) 

-1 /(« - 

l/((u-mf.)(s- m 2 h )) 
l/((u ~ m\)(s - m 2 H )) 

-1 /(s - m£) 

-l/(s - m|) 

E (flfafaxfa + fafiTuxTu + 2f[faT l xT ll ) + 2^ (fafafaxT ]U + 

i,j =1 2=1 

fafafaxfav + /i/s^xTy + JtfaTuxTm + faffafaxJ^y + ttfaT u xTy) + 
fi%\\X / T\\\ + flfayXTiy + f 2 TyXTy + 2f :i f 4 T m xT]y + 


2fafaTmXTy + 2fif 5 TiyXTy 


(All) 
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tit\ —» AA 

I. t-channel 1 2 exchange 

II. u-channel 1 2 exchange 

III. point interaction 

IV. s-channel h annihilation 

V. s-channel H annihilation 


fl 

= 

r? - 


h 

= 

r? - 


fs 

= 

C'h-h-A-A 


u 

= 

Cti-ti-h Ch-A-A 


f§ 

= 

Ch-ti-H Ch-A~. 

4 

A-x-T, 

= 



TnxTn 

= 

1 liu-mlf 


TiiixTm 

= 

1 


T w xT w 

= 

1 l(s ~ ml) 2 


T v xT v 

= 

l/(s-m 2 H ) 2 


%xTn 

= 

l/((t - ml)(u - 

m i)) 

TixTw 

= 

1 /(t - m} 2 ) 


%xTiy 

= 

— 1/ ((t ~ m 2 2 )(s 

- m 2 h )) 

%xT y 

= 


~ m 2 H )) 

TnxT m 

= 

l/{u-m 2 h ) 


AjxAv 

= 

-l/((w-m|)(s 

~ m 2 h )) 

TllXTy 

= 

-l/((t-m 2 h )(s 

~ m 2 H )) 

T m xT w 

= 

-l/(s - ml) 


T m xTv 

= 

-1 /{s~ m 2 H ) 


AyXTy 

= 

1 /(( s -m 2 h )(s- 

m 2 H )) 

T 2 

= 

fi%xTi + f 2 r u xr u + fl 


‘IV 


+ / 2 T v xT v + 
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2/i/2^ix7ii + 2/i/ 3 7ix7ih + 2/'j f{T\X%\; + 2/i/ 5 7ix7y + 
2 / 2 / 3 ^ 11 X ^111 + 2/2/ 4 7nx7iv + 2/2 /s^ixTv + 2/ 3 / 4 7nix7iv + 
2/3/5^111x7^ + 2/ 4 / 5 7f v x7^ 


Mi —» /ia [//a] 

I. t-channel t 2 exchange 

II. u-channel t 2 exchange 


/l — Cti-h-h. Cfi-h-A [Cti-h-H Qi-t2-a] 

/2 = Ctj-ia-ft Qi-ta-a [Qi-t'2-H Qi-fe-a] 

M = l/(t-m|) 2 
TnxTn = 1 /(m —mf 2 ) 2 
7^x7^! = l/((t-m|)(M-m|)) 

|T| 2 = rtT 1 xT l + tiT u xTu + 2f 1 f 2 T l xTn 


ht\ —» //+//- 

I. t-channel &n, 2 ) exchange 

II. point interaction 

III. s-channel h 

IV. s-channel H 


(A 12) 


(A 13) 


/i 

= (C il -i i - H +) 2 

/a 

= C t 1 -ti-H+-H~ 

/3 

= Ci 1 -i 1 - h Ch-H+-H- 

h 

= Ci 1 _i 1 _ H C H _ H +_ H - 

7ix7i 

= 1/((t-m?)(t-m?)) 

TfixTn 

= 1 

T m xT m 

= V(«-m 2 ) 2 
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TiyXTiy = l/(s - m 2 H ) 2 
TixTu = —1 /(t — mf.) 

TjxTm = 1 /((s-ml)(t-ml)) 

TixT w = 1 /((s-m 2 H )(t-rn- bi )) 

%\xTm = — l/(s — m 2 h ) 

TnxTiy = — l/(s — m 2 H ) 

T m xT w = 1 /{{s-ml)(s-m 2 H )) 

\T\ 2 = E /i/^x7J + flTnxTn + flT m xT m + flT w xT w + 2 E ./l/ 2 M T + 

i,j= 1 i=l 

2 E fifsTixTm + 2 E f\h%x%j + 2f 2 f 3 T u xT m + 2f 2 f 4 T u xT w + 

1=1 1=1 

2f 3 hT m xT w (A14) 

tltl -> tt 


I. t-channel gluino exchange 

II. u-channel gluino exchange 

III. t-channel X°i, 2 ,3,4) exc h an g e 

IV. u-channel X? 1)2 , 3 , 4 ) exchange 


flL 

flR 

Ji 

K, 

TlXTy 


TnxTu 


T m xT m 


f2L = f L = (qE- 


ti-t-gJ 
R ~\2 


■f‘2R = fn ~ ( C fl _ t _gj 

+ c k-t~x °)/ 2 

K-t-*>- c L t -$)l2 

(2/9) ((/E + fR) m l( s - 2m 2 ) - f L f R {Am 2 ml + 2st + 2(m| - m f 2 - t) 2 ) 

+\[JiJr (h + + m 2 - m\))/(t - m?) 2 

(2/9)((/l + fl)'m 2 g(s - 2m 2 ) - f L f R {\m 2 t m\ + 2sm + 2(m| - m 2 - w) 2 ) 

+ y/hfR (h + f R )^m t m- g {u + m 2 t - m\))/(u - m?) 2 
(—2 (4 Ji Jj Ki Kj m^o m^o s + Ji 2 ( J 2 (—— 2m 2 i m 2 + 3m/ — 4m/m^o 
+4m/my0 — Amfm^om^o + m 2 .s — m 2 s + m t myOS — m t myOS + mjpm^os 






xr x* 
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+m‘f i t + 3m 2 t — m t m^qt + m t m^qt + m^u — m 2 u + m t m^qu — m t m^qu 
— tu ) + Kj (m^ + — 3 m t + 4m t m^q + 4m t m^q — 4m t m^om^o 

—m^s + m 2 s — m t m^qs — m t m^qs + m^om^qs — m 2 t — 3 m 2 t + m t m^qt 
+m t m^qt — m^u + m t u — m t m^qu — m t m^qu + tu)) + (A',,- (— 

—+ 3 m t + 4m t m^ o — 4m t m^o — 4m t m^om^o + — m t s 

—mt'm^qs + m t m^qs + m^qm^os + m^t + 3m 2 t + m t m^qt — m t m^qt 
A-m^u — m 2 t u — m t m^qu + m t m^ou — tu) + + 2 m^m 2 — 3 m^r 

-4m t m^o — 4 m t m^q — 4m t m^qm^o — m^s + m t s + m t :m^qs + m t m^qs 
+m^qm^os — m^t — 3 m 2 t — m t m^ot — m t m^qt — m^u + m 2 u 
+m t m^qu + m t m^qu + t-u))))/((m|o — £)(—m|o + £)) 

7iyx7iv = (—2(4 Ji Jj Ki Kj m^o m^qs + J 2 (K 2 (-m-^ — 2m 2 i m 2 + 3 m* + 4m ( 3 m^ 

—4 m^rn^q — 4m 2 m^qm^q + m^s — m 2 s — m t m^qs + m t m^qs + m^qm^os 
+rn 2 _t — m 2 t — m t m^qt + m t m^ot + m^u + 3m 2 u + m t m^qu — m t m^ou — tu) 
+</,• (m^ + 2m t - i m t — 3 m t — 4m t m^q — 4 m t m^o — 4m t m^qm^q — 

+m?s + m t m^os + m t rn^os + m^om^os — rrif t + mft + m f rn^ot + m t rn^ot 

c A j Kj K% Kj Cl c A j Kj 

—m^u — 3m 2 u — m t m^qu — m t m^qu + £«)) + K 2 (J 2 (—— 2m 2 i m 2 
+3m \ — 4m^m^q + 4 m^m^q — 4m 2 m^qm^o + — m 2 s + m t m^qs 

—rnt'm^qs + m^qm^os + rri\t — m 2 t + rn t 'm^ot — 'm t m^ot + rn^u + 3m 2 u 
—m t m^qu + m t m^qu — tu) + K 2 (rrij. x + 2m 2 jn 2 — 3mf + 4mfm^o 
+4m t 3 m^o — 4m 2 m^qm^q — m^s + m 2 s — m t m^qs — mt'm^os + m^qm^os 
—m 2 x t + m 2 t — m t m^qt — m t m^qt — m^u — 3m 2 u + m t m^qu + m t m^qu 
+tu))))/((m^o - u)(-m\o + m)) 

= (-2/27)(f L f R (2(m~ i - m 2 - u)(m 2 h - m f 2 - t) - 4(s - 2m ( 2 )(m| - m 2 )) 
+\[JlJr (/l + fR)m t m~ g {t - u) + (fl + /fl)"ig(s - 2 m 2 ) - 4m 2 tf L f R 

—m t mg\Jf L f R (f L + f R ){t + u)~ 4 rn 2 m 2 gf L f R )/((t - m 2 )(w - m~)) 

Mu = 0 

TJxT^v = 0 
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0 


Tiixr m = 

T u xT w = 0 

7nix7iv = (8 J t Jj K, Kj rn^o rn^os + J 2 (—(K 2 (2mf i + Am^m 2 — 6 mf — 8mfm^o 
+8m t 3 m^o + 8m 2 m^qm^q + 2m|s + 2 m 2 s + 2m t m^qs — 2rn t rn^o s 
—2m^orriyOS — s — 2m^t — 2 m t t + 2m t m^qt + 2 m t m^ot -f t ~ 2 m^u 
—2m t u — 2m t m^qu — 2m t m^qu + u )) + Jj (2+ Am^rn^ — 6 m t 

o o o 0 0 

—8m t m^o — 8m t m^q — 8m t m^om^q + 2 m^s + 2 m t s + 2m t m^qs + 2m t m^os 
+2m^om^qs — s — 2 mf — 2 m t t + 2rn t rn^ot — 2m t m^qt + t — 2 m^u 
—2 m 2 u — 2m t m^ou + 2 m t m^ou + u 2 )) + K 2 (K 2 (2m~ i + Am^m 2 — 6 m* 

+8m t m^o + 8m t m^q — 8m t m^qm^q + 2 m^s + 2 m t s — 2m t m^qs — 2m t m^qs 

0 0 0 0 0 

+2m^qm^qs — s — 2 mf — 2 m t t — 2m t m^qt + 2 m t m^qt + t — 2 m^u 

—2 m t u + 2m t m^qu — 2m t m^qu + u ) — Jj (2+ Am^ i m t — 6 m t 

o o 0 0 0 

+8m ( m^o — 8 m t m^> + 8m t m^qm^q + 2 m^s + 2 m t s — 2m t m^qs + 2m t m^qs 

0 0 0 0 0 

—2m^qm^qs — s — 2 mf — 2 m t t — 2m t m^qt — 2 m t m^qt + t — 2 m^u 

—2 m 2 u + 2m t m^qu + 2m t m^ou + w 2 )))/((m|o — t)(m|o — -u)) 

4 

|T|- = TfxTJ + 7 Jix7Ji + (7^nx7Jii + TJyxTJy) + 27Jx7Ji + 

*J=i 

4 4 

2 ^ (7Jx7Jn + TrxTrv + ^nxTjn + 7nx7f v ) + 2 ^ 7^nx7Jv (A15) 

*=i *j'=i 

X*i —> tg 


I. s-channel t annihilation 

II. t-channel t\ exchange 


flL 

flR 

f2L 

f2R 


r L 

L 'h -t-x? 


fin 


/^iR 

G h-t-xS 


a- t - s 

Qi-ti-g 

^h-h-g 
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7M = (-4/6)((/ 1 2 l + /iij)2((3m 2 - s)(s + m\ - m\) + (s - m 2 )(m| + m 2 t - t)) 

—8/il/ir m t mjj(i + m 2 ))/(s - m 2 ) 2 

TuxT n = (-4/6)(2 (t + ^| 1 )((/|l + fm){ m \ + "i? ~ *) ~ 4m^mt/ 2L / 2R ))/(t - mfj 2 
= (—4/6 )((/i L / 2 l + /iR/ 2 i?)/ 2 ((m| + m 2 - t)(s - 2m| + 2m| + m 2 ) 

— (s + m 2 )(2s — 3m| — 2m 2 x + w) + (2 ml + 3 m 2 — 2-u — s)(s + m| — m|) 
+2m 2 (2s - 3m| - 2m| + «)) - (/il/ 2 r + /ir/ 2 l) m t m^(4m| + 4m 2 
-2m| -2u))/((s- m 2 )(t - mfj) 

|T| J = 7Jx7J + TjjxTJj + 2TixTu (A16) 

Xti —» tz 

I. u-channel 234-) exchange 


Pll 

Plr 

Prl 

PrR 

PxT, 


T 


2 


C L 

G h-i-x° 

r^R 

G h-*-x° 

riR 

G h-*-x° 


C|o .-.0 ^ 

Xi-Xi“ z 

r^iR 

Xi 

C ~0 ^0 7 

Xi-Xi“ z 
siR 

O~o ^0 7 

Xi“Xi-X 


(!/ 2 )((/ll /ll + /m?. Prr)H u - m h + )(« + m x “ m x) “ ( m x + - *) M ) 


+((m| + m| — w)/m|)((m| — m| — tt)(w — m| + m 2 ) — (s — m 2 — m 2 z )u)) 
HPlr Prr + Prl /ll)wtc*°((« + m| - m|) + (l/m|)((m| + m| - u) 
x (m| - m| - «))) + (P LL f RL + /rr p LR )m t m^o((u + m\ - m 2 z ) 

+(l/m|)((m| + m| - w)(m| - m| - u))) 

HPlr Plr + Prl /rl)™x°™x°((* + m ! “ m t) 

+( 1 A n l)((™* + m| - «)(s - m 2 - m|))) 

+ (/ll Prr + Prr flLprutm^u + (f LL f LR + f l RR P RL )3m-o(u + m 2 - mfjm* 

+ (/lr /ll + Prl Prr)S m^u + m 2 - mfjm* 

+ (/lr Prl + Prl ~ ~ m I°)) 

E ( Al7 ) 

*J=i 
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xti —> bW+ 

I. u-channel jP 2 ) exchange 


ri _ /^iL riL 

III ~ °t 1 -6-x+ L 'x° 1 -x~-w+ 

fi _ /~<L /~iR 

Jlr - ^ii-b-xt Ly x° 1 -x~-w+ 

fi _ s~iR riL 

Jrl - ° h-b-xt Ly x° 1 -x~-w+ 

ri _ s^iR f'lR 

Jrr - Ly x° 1 -xi-w+ 

T i xT i = (1/2 )({P LL f LL + f RR f RR )(((u - m\ + mg) (u + m| - m 2 w ) - (m| + mg - t)u) 
+ ((m| + m 2 w - w)/m^)((m| - rn 2 w - u){u - m| + m 2 b ) - (s - rn b - m 2 w )u)) 
+(Plr Prr + Prl /i L )m 6 m-+((u + m| - m^) 

+ (l/m^)((m| + - u){m\ - rn 2 w - u))) 

HPll f 3 RL + Prr /i R )m 6 m^+ ((m + m| - m^) 

+(l/m 2 w )((m\ + - u)(m| - m^ - «))) 

+(/lr Plr + Prl /^ L )m-+m-+((t + m| - mg) 

+(l/m^)((m| + m 2 w - u)(s - m 2 b - m 2 w ))) 

+ (/ll Prr + Prr pLLp^bm^u + {P LL f LR + P RR f RL )3m~+(u + ml - mfjm* 


+(/lr /ll + Prl P R Rpm~+( u + mg - mfjm* 

HPlr Prl + Prl /i^6m~+m~vm fe m*)/((t - m|+)(t - m|+)) 

EM 

*j=i 


(A18) 


—- tH [th] 

I. s-channel t annihilation 

II. t-channel t(i t 2 ) exchange 

III. u-channel X°i, 2,3,4) exchange 


A L = CX- S C t - t -H CU-kl 

Ah = Cf_,_ s C,-,- H [c£_,^, C,.,. h ] 
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f l 2L 

fm 

f 3 LL 

f:SLR 

I'Irl 

f'iRR 

%xT Y 

T u xT u 

T m xT m 


TixT u 


%xT m 


TnxTm 


C h-t-x\ C h-h-H 

CLt-t ? C h-h-H 

pyL 

G il-t-X° G X?-X°-^ 
s~iR syL 

G il-t-x° G x?-X°-^ 

a'y/? 

G ii-t-x° G x?-x°-h 


C h-ii-h. 


[Cg+x ? 

Pt*-x? C S-x?-fc] 

[^-t-x? 

[Cg-^-x? ^-X?-J 

(V 2 )((/il + /ir)((s + "4 - m 2 H[h] )(s + m| - m|) - 


-t) + 2 m 2 t (s + m| - m|)) - f 1 L fiR^m^m t ( 2 s + 2m 2 


( 2\ / 2 i 2 

(s — m t ) (m£ + m t 


O-ftWLfL + f2RfL)( m2 x + m t t) 2m t m x(/2i/2R +/m/Il)) 


(1/2)((/3Z,k/3LR + f3RLf3RL)(( m t m fi+ M )( m | + M ) 

+ m 2 - t)) + {f 3LL f£ LL + r 3 R R fmR}rn^m^(ml + m 2 - t) 
~(f 3 L R f3LL + f 3R Lfi RR ) m x°( m t _ m fi + M ) ~ (/ 3 LL/ 3 ZJ? + fl RR fi R L) 

x m^o m*(m 2i - m\ + u) + {f 3RL f J 3LL + f 3LR f J 3RR ) - m| [h] 

+ w ) + ( P 3RR flLR + fl llHrl) m t m x°S m \ - m H[h] + u )~ (/m/L 
+ /3Li?/sRL) 2 ™X M - UIrrILl + PllMrr) 2 ™x° m x° m x) 

4 J 

/((« - m |o)( M - m l°)) 

(1/2)((/il/2l + fiRf 2 R )m t (s + 2m| + m 2 - m| - t) - {fi R f l 2L 
+fiLf2 R )mx{s + 3m 2 - m| w ))/((s - m 2 )(t - m|)) 

(l/2)((f 1L f 3RL + /ifl/3Lfl)((s + - m|)(M + m 2 - m?) - (m| + m 2 

-t)(m 2 + m| - m| - m^ [h] ) + (w + m| - m^ [fc] )(s + m 2 - m^ [h] ) 

+2mj(« + "4 - m^ [h] )) + (Al/sll + .UrPrr) ™x° m t {s + 2m| + m 2 - m| 
“4 - {flLfzLR + /ik/sRl) ™X m < ( M + 2 ™t + - 2m l ~ m H[h ]) 

~{fiLf l 3RR + fmflLL) m x m x° (S + 3m 2 - m^ w ))/((s - m 2 )(w - m 2 o)) 

(l/ 2 )((/2i?/3LR + / 2 L/ 3 RL) m < ( M + m | — m !ff[h]) + (/ 2 L/ 3 LL + fi R f 3 R R ) 

x m^o (m| + m 2 - t) - ( f 2L fk R + fmfluL ) ™x ( M + - m |) 

-(fmfsLL + fLfsRR) 2 m x m t m x°)/(( t - m |)( n “ m l°)^ 
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2 4 2 4 

\T\- = 7)x7) + 22 TnxTn + 22 ^m x ^m + 2 y~)T^x^n + 2 22 A x ^m + 

*J=1 i,j=l 1=1 i=l 

2^^ 7fix7rn 

j=i *=i 


(A19) 


xh —> tA 

I. s-channel t annihilation 

II. t-channel t 2 exchange 

III. u-channel X(*i 2 3 4 ) exchange 


c Lt-*° c t - t -A 

c L t - x ? 

C h-t-x° C h-t2-A 
C Lt-x\ C h-h-A 

s~i L siL 

G ti-t-x? °X?-X°-A 

/oL 

G x?-x°-a 
y^fA y^fL 

G h-t-Xi ^x\-X°i-A 


7hx7h = 


fi _ /~iR /~iR 

hRR ~ Wi-t-x° ^xO-x^-A 

7^x71 = (1/2)((/^ + A 2 r ) ((s - m 2 - m^)(s + m\ - m|) - (s - m 2 ) (m| + m 2 - t)) 
-fiLfiR 4 m* m t m 2 A )/(s - m 2 ) 2 

7)ix7)i = ( 1 / 2 )((/ 2 l/ 2 L + / 2 r/ 2 r)(i«| + m 2 — t) — 2 m t ^(AlAr + A?Al)) 

/((*-<) 2 ) 

AiixTiii = (1/2 )((P 3LR f 3LR + PsrlPsrl) ((™t ~ ml + «)(m| - + «) - «(m| + m 2 

+ (PllMll + IIrrSIrr)™x° m^m\ + m 2 - t) - {P 3LR P 3LL 
+/Al/AA™x m x °( m t - m l +u)~ (P 3LL fi LR + PshhUrl)™^ ™x( m ? 
-m 2 h + u) + (. P 3RL f 3LL + f 3 LRf3RR,) m t m x° ( m i - m A + «) + ( PrrJLr 
+P iTjJlRL)m t /m^im\ - + w) - (P 3RL f 3LR + PlrPrl) 2 m t "A M 


-( PrrMll + PsLLfsRR ) 2 m^o m^o m t m*)/((u - m| ? )(« - m 2 o)) 


7tx7tt — 


( 1 / 2 ) ((/WA - f2LpL)‘m t {u - m\ - m^) + ( fnf\R ~ Aa/il) ™x( s “ m t 
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-m\))/{{s-m 2 t ){t-m h 2)) 

TixT m = (l/2)((f 3LR fiR - IIrlUl) ((m| + m t 2 - t) (rn 2 + m| - m| - m^) 

— (s — rn 2 — m 2 A ){u + m| — m^) — (w + rn 2 — m~ i )(s + m| — m| i ))/2 
HflR R fiR ~ fluJuf) m^o m t (u - m\ - m|) + (f 3LR fiL ~ /Wik) 
x m t m^(u - m\ + m^) + (. fl RR fiL ~ fl L LhR)rn^mhm 2 t + - s)) 

/((s-m 2 )(w-m|o)) 

^n x ^m = (l/ 2 )((f 2 R f- 3L R + /jl/Zrl) m t (u + m | - m^) + (/2L/3LL + ImflRn) 

x m^o(m| + m 2 - t) - (/ 2 l/slr + f2Rfl RL ) m h u + 

HhRfsLL + f2Lf3RR)2m^m t m^)/((t - m\){u - m 2 o)) 

4 4 4 

|T|~ = 7ix7i + TnxTn + ^ 7fiix7f IT + 27fx7f T + 2^Ax^ni + 2 TrixTyn (A20) 

i,j=l i =1 i=l 

Xh —» bH+ 

I. s-channel t annihilation 

II. t-channel x A 2 ) exchange 

III. u-channel &( li2 ) exchange 


/ill 

nL 

(-iL 

°ti-*-x? 

flLR 

- c L 

— L t-b-H+ 

riR 

°h-t-x? 

flRL 

s~iR 

(~iL 

°h-t-x? 

flRR 

/—iR 

/—iR 

°tl-t-x? 

/ 2 LL 

= C; h -+ 

ti-b-xj 

riL 

Ly x° 1 -xt~R + 

/ 2 LR 

= C; h -+ 

ti-b-xt 

iR 

x?-x+-R+ 

f'2RL 

/—iR 

- W 7 , — 1 - 

ti-b-xj 

C L 

°x?-x/-R+ 

fmR 

(-iR 

- V_X 7 , - 1 - 

ti-b-xj 

&& 

1 

Xi 

*•+ 

1 

fc! 

+ 

fk 

- C *- L 

- 

C il _- b ._H+ 

fH 

/-iR 

- °Si- 6 -x? 

C il _~ bi _H+ 


^i x ^i — ( 1 / 2 )((/il_r + /i.rl)(( s + 


m R+)( s + m \ ~ m \) - s(m| + m 2 - u )) 
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+ {fiLL + fiR R ) m t( m l + m b ~u)~ 2(f 1LL f 1LR + fiR R fiR L ) m^mt 
x 0 + ml - m 2 H +) + 2(f 1LR f 1RR + fiLLfiRL)m b m t (s + m| - mfj 
~4(/i LR /i RL )m & s - 4:f 1LL f 1RR m 2 t m* m b )/{s - m 2 t ) 2 
T n xTn = (1/2)((f 2flJ ^ RR + /sll/sllX^ + - ™|)(* + - ™h+) - + ml 

“«)) + (f l 2RLf-2RL + ILrILr) \+ ( m l + “ M ) “ (KrrArL 

+f2LLf2LR) ™X m x+ + m b ~ m h ) “ ( P 2 RLP 2 RR + ILrILl) m * 

x (t - ml - m\) + {f\ LL f 2RL + f l 2RRf 3 2LR ) ™ b m^+ (t + m\ - m 2 H+ ) 

+ (/2Z,il/2i?R + /Irl/IlZ,) m b m x+ X + ^ ~~ m H+) ~ (f^LL^RR + f^RR^LL,) 

x 2 m b rri^ t - {f 2LR f 2RL + /a rlUlr) 2 m x+ m x+) 

/((t - m\+)(t - m 2 + )) 

A. x A.j 

= ( 1 / 2 ) ((/ 3 Z./ 3 Z, + /3R/3ii)( m l + rn t — u) — 2 m^m b {fl L fi R + fl R fl L )) 
/{{u-ml){u-ml)) 

TixTn = ( 1 / 2 ) {(r 2RR fi LR + r 2LL fiRL)(l/mt + ml - m\)(s + m| - m?) - (m| 
+m£ - w)(m| + m 2 b - - m#+) + (s + m/ - m^+)(t + m| - m#+)) 

HfwLfiLL + f 2LR fiR R ) m~+ rn t (m| + - u) - (f 2RR fiLL + P 2 uJ\nn) 

x m t (t + ml - mfj - (f 2RL fi LR + f 2LR fiRL) m^+ m * (s + - m|+) 

+ (/2Ll/iRR + f2RRflRR) m bm t (t + J71- — m 2 H+ ) + (f 2 L R flLR + /Irl/i-Rl) 
x m b m^+ (s + m| - m|) - {f 2LL fiLR + f 2RR f\RL) m b m .* (m| + m 2 b - m\ 
~m 2 H + ) ~ if 2LR fiLL + r 2RL fiR R ) 2m b m^m t m-+)/((s - m 2 t )(t - m|+)) 
^ix7J n = (l/ 2 )((/iLfl/ 3 fl + /irl/ 3 l) m b (s + m| - m|) + (/ill/^ + fi RR fl R ) m t 
x (m| + m 2 b -u)- (/1 lrPzl + I\RtJIr) m% {s + m/ - m^+) - (Jill fan 
+fiR R fl L ) 2 m b m^)/((u - m 2 bi )(s - m 2 t )) 

TiixTiii = (1/2 ){{f 2LR fl R + /Irl/sl) m b (t + m| — m 2 H +) + (/ 2 LL/ 3 L + XIrrIIr) 
x m x+( m l + - M ) - (/Ilr/sl + XIrlIIr) mx (t + m 2 b - m|) 

-(/ 2 LL/ 3 R + fuwfk) 2 m x + m x)/(( M - - m l +)) 

L A.J 

2 2 2 2 

\'T\- = 7Jx7I + ^ TrixTii + ^ ^iiixTJh + 2^TJxTJi + 2^7Jx7Jn + 
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(A21) 


2 

2 JZ ^n x ^ni 

i,j= 1 

—■» t£ 

I. t-channel 234) exchange 


ri _ a~yL 

Jll - °h-t-x° Wh-^-x? 

ri _ /^L 

Jlr - °ti-t-x° 

JRL - H-t-x? °€i-y-x? 

Jrr - H-t-x? 

71x71 = (( P LR f RL + /rl/lrXC^I - "X “*)(* + "»? - "ifj - t(s - m\ - m})) 

+ (/1r/lr + /rr/rr) ™x? m x° ( s “ m t “ m T) 

-(/lr/ll + /rl/rr) ™A™x° ( m fi - m t 2 - *) 

-(/ll/rl + f RR f L R ) rn t rn^ ( m ? - m 2 - y) 

+ (/rl/ll + /lr/rr) "»t "1*0 0 + "1? - "if,) 

+ (/rr/rl + /ll/lr) "it (* + "if - "if,) 

-(/rr/lr + /lr/lr) "it m t 2 1 - (f RR f LL + f LL f RR ) 2 m*m £ m*o m*o) 

/((t-m|o)(t-m|o)) 

|T| 2 = (A22) 

*j=i 

t\I —■> 6 z/ 

I. t-channel X(i 2 ) exchange 


/ll 

= -+ 
ti-6-Xf 

M * t"* 

1 

1 

•*• + 

Plr 

= C'f -+ 

ti-b-xj 

1 

1 

-• + 

/rl 

r^iR 

— w r » ~ + 

tl-b-xj 

,0 

1 

1 

-• + 

/rr 

(~iR 

°h-6-x+ 

,0 

1 

-•+ 
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M = (( P LR f RL + /rl/lr) ((ra| - - t )( t - m \)- t ( s - m 2 b )) 

+ (/1l/lL + /rr/rr) m x + m x + ( s — m b ) 

+(/rl/ll + /lr/rr) ™b m ~+ (t - m \) 

+ (/rr/rl + /ll/lr) ra fe (* _ m !i)) 

|T| 2 = EM (A23) 

i,i=l 

tl/* —■> 

I. t-channel X(i, 2 , 3 , 4 ) exchange 


/ll 

/lr 

f l RL 

/rr 
%x 71 


T 


2 


/'-(R 

0 ii-t-X° °h-^-X° 

/OR /^L 

/OR S~lR, 


((/1l/lL + fRRfRR)(( m l ~ “*)(* + - "1?,) - *0 - m t ~ m2 ()) 


+(/k/k 

+ f l RL 

/lr 

) m x 

r m x 

o (5 - 
/ 

- - 

- ra 2 

— (/ll/rl 

+ /rr 

1 /lr 

:)ra* 

m { o 
Xj 

(ra| 

-ra 2 

-*) 

— (/lr/ll 

+ /rl 

/rr 

) 

m { o 

(ra| 

-ra 2 

-*) 

+ (/rr/rl 

+ /ll 

/lr 

)ra t 

m,-o 

Xj 

(t + 

ra 2 - 

ra|) 

+ (/rl/1l 

+ /1r 

/rr 

) rn t 

m^o 

A-i 

(t + 

ra 2 - 

ra|) 


(IrrIll + /ll/rr) m t m e 2 1 (/rl/rl + /lr/lr) 2 ra t ma m^o o, 

/((t-m|o)(t-m| 0 )) 


E M 

*J=i 


(A24) 
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